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Abstract. I study the Lyapunov exponent and the integrated density of 
states for general Jacobi operators. The main result is that questions about 
these, can be reduced to questions about ergodic Jacobi operators. Then, I ap- 
ply this to a{n) = 1 and b(n) = /{nf (mod 1)) for p > not an integer, and to 
obtain a probabilistic version of the Denisov— Rakhmanov-Remling Theorem. 



1. Introduction 

This paper is part of my effort to study the Schrodinger operator, 

(1.1) {Hu){n) ^ u{n+l) + u{n- 1) + /{nP (mod l))u(n), 

where u(—l) = 0, / : [0, 1] ^ M is a continuous function, and p > is not an 
integer. I was intrigued by the fact that for < p < 1 and /(O) ^ /(I), one 
has the absence of absolutely continuous spectrum and vanishing of the Lyapunov 
exponent on an interval. This is somewhat surprising since (mod 1) has nice 
uniform distribution properties. We will discuss properties of these operators in 
Section [5l In particular, we resolve the discrepancy between the perturbative and 
numerical calculations of Griniasty and Fishman in [6J in Corollarv l5.2l bv proving 
an exact formula. 

In order to understand the consequences and reasons for zero Lyapunov expo- 
nent, it turned out to be useful to work with general Jacobi operators, which are 
introduced by 

J : -^£^m 

(1.2) 

Ju{n) — a{n)u(n + 1) + h{n)u{n) + a(n — l)u(n — 1), 

where Cq"^ < a{n) < Co and —Co < b{n) < Cq for some Co > 1- We let m±{z) 
be the Weyl-Titchmarsh m functions of the restrictions of J to ^^(Z±). J is called 
reflectionless on A if 



(1.3) m+{t) = -TO_(t) 

for almost every t & A. Denote by L{E) the Lyapunov exponent of J, by J^") the 
n-th translate of J, and by 6j the Dirac measure. We have that 

Theorem 1.1. Assume L(E) = for almost every E ^ A and 

(1.4) = ^i-^ ]^ E '5,/<". 

n=0 
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in the weak * topology, then /i almost every Jacobi operator has absolutely contin- 
uous spectrum in the essential closure of A and is reflectionless there. 

Remling has shown in [12] a very similar result. He has assumed that A C a^dJ), 
and concluded that every J in the uj limit set of J*^"-* is reflectionless on A. Since 

(1.5) (Tac(J)C{S: LiE) = 0}, 

the assumptions of the above theorem are weaker, but also the conclusion is. One 
can easily check that sparse potentials (as discussed in [T2]) provide examples, that 
show that this distinction is sharp. 

The above theorem will follow Theorem 14.11 which provides a formula for the 
Lyapunov exponent L{E) in terms of the Lyapunov exponents of the ergodic families 
arising in the ergodic decomposition of the limit measure /i. 

Theorem 11.11 implies in particular the following result, which has to be thought 
of as a probabilistic analog of the Densiov-Rakhmanov-Remling theorem 

Theorem 1.2. Let J be a Jacobi matrix with acss{J) — [^2,2] and L{E) — for 
almost every E e [—2, 2], then for every e > 

(1.6) lim 4#{1 < 71 < TV : \a{n) - l\ > e or \b[n)\ > e} ^ Q . 

We will obtain a generalization of this theorem for finite gap operators (see 
Corollary I4.3|l and the above claim for the Lyapunov exponent as corollaries of 
Theorem 14. II in the next section. 

In Section[21 1 collect a few results on the space of all Jacobi operators and discuss 
measures on that space. In Section [3l I discuss results about ergodic Schrodinger 
operators. The main results are stated in SectionlH The application to the potential 
V{n) — f{nP (mod 1)) is examined in Section [5l Section [6] proofs some facts 
about the Lyapunov exponent for general Jacobi matrices and provides the proof 
of Theorem 14. 1[ which has to be considered the main result of this paper. 

2. Probabilistic averages of Jacobi matrices 

Given bounded sequences a : Z ^ (0, oo), 6 : Z — > M, we introduce the associated 
Jacobi operator J : P{Z) —> i?^(Z) by 

(2.1) {Ju){n) = a{n)u{n + 1) + b{n)u(n) + a{n - l)u{n - 1). 

We will often identify J with (a, b). Fix now Co > 1, and introduce as the set of 
all Jacobi operators, such that (a, b) satisfy the inequalities 

(2.2) < a{n) < Co, -Co < b{n) < Co- 

We endow with the strong operator topology, which just corresponds to pointwise 
convergence on the level of the sequences (a, b). We remark that J' is now a compact 
metric space, where an explicit example of the metric is 

•^~) = E (J - J)Sn)\ + \{Sn, (J - J)Sn+l)\) 

riGZ 

(2.3) = J2 - + l«W - 
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Denote by S the shift operator on £^(Z) that is 

(2.4) {Su){n) =u{n + l). 
Introduce S : J ^ J by 

(2.5) SJ = S*JS, 

and denote = S"' J. 

We will denote by the space of all Borel probability measures on J. For a 
Jacobi matrix J, we introduce the corresponding Dirac measure Sj by 

\ J e A 



(2.6) 5j{A) = 



J iA. 



We will be mainly interested in the limit points of the averages of the Dirac measures 
of the translates of a Jacobi matrix. For this, introduce for a Jacobi matrix J and 
an integer > 1 the average 

(2.7) A^,j^^Y.5j,.,, 

71=0 

which will be a measure in A^^. 

We denote by lli{J) the (topological) ut limit set of the translates of J, that is 

(2.8) uj{J)^{JeJ: 3nj^oo: J = lim J*"^'}. 

We recall that the weak * topology on gives rise to the notion of convergence 
fJ-n —>- fJ- for every continuous function f : J' we have 

(2.9) lim / f{J)iin{J) = / f{J)KJ)- 
We remark 

Lemma 2.1. A4^ is a compact and metrizable space in the weak * topology. 

We write supp(/i) for the support of a measure /x, which is the smallest closed 
set A, such that fJ.{A) = 1. Furthermore, we call a measure n G Ai^ shift invariant, 
if for any Borel set A C J' 

(2.10) fi{A) = fi{SA). 

Lemma 2.2. If fJ. = limj^oo A^.^j in the weak * topology, then 

(2.11) supp(m) Ccc;(J), 
and fj, is shift invariant. 

This implies the following consequence 

Lemma 2.3. Let fi = limj^oo A^. (J) for some Nj — > oo. Then for fi almost every 
J , we have that 

(2.12) (7ess(J) C CTcss(J). 

Proof. Follows from the fact, that the inclusion holds for every J G lu{J). □ 

We will need the following result from measure theory, known as Portmanteau- 
Theorem (see e.g. [SJ Theorem VIII.4.10.]). 
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Theorem 2.4. ^ /i m the weak * topology, is equivalent to that for every Borel 
set B with fJ-{dB) ~ 0, we have that 

(2.13) lim Hn{B) = fi. 

n — >oo 

We say that a sequence J*^"^ converges to a set A C along Nj oo in 
probability if for every e > 



(2.14) Hm -^#{1 < n < iVj : d{J^"\ A) > e} ^ 0. 

We have the following resuh. 
Lemma 2.5. Let 

(2.15) ^i = lim An^j, 

and S a support for fi. Then J'^"-' converge to S along Nj in probability. 

Proof. For e > apply the last theorem to i? = { J : dist(S', J) > e} to conclude 
that ANj,j{B) ^ as j — > cxD. By rewriting, one sees that this is exactly the 
definition of convergence in probability. □ 

For A C Z, denote by J a the restriction of J to ^^(A). We will caU f : J ^ R 
compactly supported, if there is a finite set A C Z such that 

(2.16) f{J)^f{h 
whenever Ja = Ja- We have that 

Lemma 2.6. For /i„,/i G , we have 

(2.17) fin M 

in the weak * topology, if and only if for every compactly supported f 

(2.18) lim f fdfia = f fdfi. 



Proof. It clearly suffices to show that p.lSp implies weak * convergence. So given 
/ : — > M and e > 0, we have that show that there is an iV > 1 



yn> N 



J fdfJ-n - J fdfi] < e. 



Since, / is continuous, we may find for each J E J an integer K > \ such that 
|/(J)-/(J)| <|,where 

J e UjJi = {J : J[-K.K] = J[-K,K]}- 

Since the Uj^k are open sets and J7 is compact, finitely many of them cover J7. In 
particular, we can find a maximal necessary K . Hence, we may approximate / by 
/, which is supported on [—K,K] such that ||/ — /||oo < §■ Now the claim follows 
by ([218]) . □ 
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3. Families of ergodic Schrodinger operators 



In this section, we collect basic facts about ergodic Jacobi operators. For the 
Jacobi operator background see [3] or Section 7 of 13J. For the measure theoretic 
part, see [7] or [9] 

Denote by M.g the set of all shift invariant measures. One can check that A4g 
will be a convex set, and in particular we will write £ for its extremal points. 

It is a known fact in ergodic theory, that £ are exactly the ergodic measures of 
the dynamical system 5), where 



(3.1) 



S{J) = S*JS. 



Furthermore, it follows from Choquet's theorem that one can write any measure 
/i S as a generalized convex combination. That is, there exists a measure a on 
£ such that for any f : J ^ M. continuous, one has 



(3.2) 



fdp da{(3). 



We call an ergodic measure /? G 5 on the space of Jacobi operators J a family 
of ergodic Jacobi operators. One knows from general fact, that there is a set S(/3) 
such that 

(3.3) S(/3) = a(J) 

for [3 almost every J. We may define its Lyapunov exponent by 



(3.4) 



7/3 (^) 



lim — 



log 



J 



n 

= N- 



1 



z - h{n) 1 
^ a{n) V a{nf 



d/3(J), 



and its integrated density of states by 



(3.5) 



kp{E) 



J 



(^o,X(- 



(J)5o)dp{J). 



We note that this quantity is equal to 

1 



(3.6) 



kp{E) - lim - / ir{P(_^^E){J[o,N-i\))dl3{J). 

IS -^oo 1\ J j- 



We will need the following result of Kotani theory. 
Theorem 3.1. Denote by Z the essential closure of the set 
(3.7) {EeR: jf3{E) = 0} 

then [3 almost every J has purely absolutely continuous spectrum on Z and it is 
reflectionless there. 



We define 
(3.8) 



log(A^) = / log(a(0))d/3(J). 



J 



Furthermore m^{z,J) = ((5o, ( J+ — z) ^Sq), where J+ is the restriction of J to 
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Lemma 3.2. We have that 
(3.9) 



7/3(^) = \og{A^') - / log|m+(z, J)|d/3(J) 
Jv 

= log{A/)+ / log \t - z\diy,3{z) 



(3.10) 

for every z G C. 

(|3.10p is known as the Thouless formula. It implies that 



(3.11) 



lim-fpiE + ie) = 7^(£') 



for every i? e R by monotone convergence. 

We now make the connection, to the usual definition of ergodic Jacobi operators 
(see also Section 2 in T). Let (17, T,^) be an ergodic dynamical system, and 
a : rt —> (0, oo) and 6 : f2 — > R are measurable functions satisfying 

(3.12) < a{uj) < Co, -Co < b{uj) < Co 
Co 

for almost every w. Then we can define a map 

(3.13) f-.n^j 

by f{oj) being the Jacobi operator with coefficients 

(3.14) {(a(r"c^),5(T"c^)}„ez. 
Introduce a measure (3 oti J given by 

(3.15) p{A) ^ p.{f-\A)) 

for Borel subsets A Q . Then the usual definitions of the Lyapunov exponent and 
the integrated density of states will just be 7^ and kp. 

4. Statement of the results 

For a Jacobi operator J ^ J and a sequence Nj — > (X), we introduce the Lyapunov 
exponent by 



(4.1) L(z, J, {N,}) = limsup log 



n 



1 



z — b(n) 1 



where J = (a, 6), and the integrated density of states by 

(4.2) k{E,J,{N,}) = hm ±triP^_^EMo,N,-i])), 

where the limit is assumed to exist. 

Theorem 4.1. Assume that we have a measure a on £ such that for Nj 
(4.3) 



hm -L ^ Sji„y = I pdaiP) 



then 
(4.4) 



L{z,J,{N,}) 



-fp{z)da{f3) 
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for Ini(z) > and almost every z G R, i 
(4.5) kiE,J,{N,})^ I kpiE)da{P). 



The statement is actually stronger, since one may replace the almost every by 
quasi-every (in the sense of potential theory). Theorem ll.ll is now an easy corollary. 

Corollary 4.2. Let A C R &e a set of positive measure. Assume that 

(4.6) L{E,J,{N,})^0 

for almost every E £ A. Then fi almost every J is refiectionless on A. 

Proof. Since "(p{E) > for every (3, we can conclude by (|4.4p that for a almost 
every /3, we have 

ME) = 

for almost every E E A. The result now follows by Theorem 13. II □ 

We also obtain the following probabilistic version of the Denisov-Rakhmanov- 
Remling theorem (see [1], [IT], [H]). For this recall, that a set e is called a finite 
gap set, if 

(4.7) e = [Eq,Ei] U[E2,E3]U---U [E2g+i, E2g+2]. 

Furthermore, one has a finite dimensional torus T(e) of refiectionless Jacobi-operators 
which have spectrum e. We have that 

Corollary 4.3. Let e be a finite gap set, and assume that 

(4.8) ae..(J) = e 
and L{E, {Nj}, J) = for almost every i? G e, then 

(4.9) J(") -> r(e) 

in probability along Nj, where T(e) denotes the isospectral torus. 

Proof. Assume there is a subsequence of Nj such that convergence in probability 
does not hold. By passing to a further subsequence of Nj, we may assume that 
(j4.3p holds. Since 7/5(2;) > everywhere, it follows from our result that for a almost 
every /3, we have that 

S(/3) = e 

and that 7^3 = on e. Hence, Kotani's theory implies that these operators are 
refiectionless on e, which implies in turn that /3 is supported on T(e) (see e.g. 
Section 8 in [U]). This is a contradiction by Lemma [2.51 □ 

5. The family of potentials 

In this section, we examine the family of Schrodinger operators given by (jl.ip in 
some detail. Introduce for a continuous function / : [0, 1] ^ M and r < p < r + 1, 
where r is a nonnegative integer, the sequences 

(5.1) a{n) = 1, b{n) = f{nP (mod 1)). 

Denote by J the associated Jacobi operator. For a G [0, 1]\(Q, introduce the skew- 
shift T„ : [0,1]'' ^ [0,1]'' by 

(5.2) (r„w)fe = < 1 / I. / 1 

LOk + Wfc_i 1 < fc < r - 1. 
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Similarly as in the last part of Section [3] we let (3a be the measure on J given by 
the pushforward of the Lebesgue measure on [0, l]*" under 

(5.3) [0, lY3u^ {1, f{{T^u)r)}n^z e J- 
Lemma 5.1. We have that 

(5.4) 1™.,1^E'^^<"' = / 



n=0 







in the weak * topology. 

Proof. By Lemma 12.61 it suffices to check convergence for compactly supported 
functions g. Since, every such g will be a continuous function of {b{''T')}n=-K 
some K > 1 it suffices to check that these converge. Since, / is also continuous. It 
suffices to show that 

has the same distribution in [0, 1)2^^+1 ag the orbits of the skew-shifts would as 
n ^ oo. 

Furthermore, by Lemma 2.1 in [8J and an easy argument we see that both {n+j)P 
and (T"'^^u!)r are essentially given by degree r polynomials, and thus uniquely 
determined by 

{(n + jr}-=o and {{T-+^Lu)rV^^,. 
Now Lemma 2.3 in [5] implies that the coefficients of the first polynomial are uni- 
formly distributed, and a quick computation shows the same for the skew-shift, 
finishing the proof. □ 



This lemma combined with Theorem 14.11 implies 
Corollary 5.2. For almost every E, we have that 

(5.5) L{E) = f j0AE)da. 

Jo 

This corollary resolves the discrepancy between the numerical and perturbation 
theoretical computations in [6] and shows in particular that the Lyapunov exponent 
only depends on the integer part of p. 

In particular, in the case of r = 0, that is < p < 1, one can compute that 
7/3^ {E) = for exactly 

(5.6) i?e [-2 + /(a), 2 + /(«)]. 
Hence, we see that L{E) — for 

(5.7) [-2 + max(/),2 + min(/)], 

which was first observed by Simon and Zhu in [14 for continuum Schrodinger 
operators. We observe further spectral properties in the following result. 

Theorem 5.3. We have that 

(i) Stolz: If f extends to a smooth function on the circle, then J has purely 
absolutely continuous spectrum in [—2 + max(/), 2 + min(/)]. 

(ii) Remling's Oracle: ///(O) ^ /(I); then the absolutely continuous spec- 
trum of J is empty. 

Proof. Part (i) is [TS]. Part (ii) follows from Remling's Oracle Theorem, which is 
found in fTH. □ 
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6. The integrated density of states and the Lyapunov exponent 



In this section, we will proof Theorem 14.11 For this we have to discuss some 
further properties of the Lyapunov exponent. 

We will now assume that we are given a fixed Jacobi operator J : £^(Z) 
£2(Z). We denote by J+ its restriction to £^(Z+), Z+ = {0,1,2,3,...}, and by 
Ja its restriction to ^^(A) for A C Z an interval. Denote by Ej{A) an increasing 
enumeration of the eigenvalues of Ja, introduce the density of states measure ^'„ of 
Jlo,n-i] by 



n — 1 

(6.1) „^ = -ys 



n-1 

'iS,-([0,n-l])- 



n 

For Im(2;) > 0, introduce the Weyl-Titchmarsh m function 

(6.2) m+{z,J)^{SoAJ+~^r'So). 

We will show the following theorem, which is essential in the proof of Theorem l4.1l 
Similar results can be found in Poltoratski-Remling [lOj . 

Theorem 6.1. Given a sequence Nj —>■ oo, assume that 

(6.3) lim — V (5j(„) = fi. 

n— 

in the weak * topology on . Then 

(i) The integrated density of states measures converge in the weak * topology 

(6.4) lim vn. — v. 

(ii) For Im(z) > and A as defined in (j6.7p . we have that 

(6.5) L(z,{iV,})=log(A-i)- / log|m+(z, J)|d/i(J). 

J J 

(iii) For almost every _E G M, we have that 

(6.6) T(E, {Nj}) = YiTCL~L(E + ie, {N^}). 

Since the map (a, h) i-^ log(a(0)) is continuous, we see that ()6.3|) implies 

(6.7) A exp l^j ^ log(a(0))dAi( J)^ = /i™ E log(a(n)) 

where A is the constant in (|6.5p . We are now ready for 
Proof of the Theorem \4-l\ We first observe that 

log(A-i)=^log(^^i)da(/3). 
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We may compute for almost every E, that 
L{E) = lim L{E + ie) by ([61 



\im(^\og{A-^)-J log\m+{E + te,J)\dfi{J)^ bv (1^ 

lim (^J log(A-i) - log \m+{E + ie, J)\d(3{J)^ daifJ)^ by (gj) 
lim f / 7^(i; + ie)da{(3)] by (EH) 



7;3(i?)da(/3) by(E3Il). 
This implies the first claim. The second claim follows by Thouless' formula. □ 

We now proceed to prove Theorem 16. II Introduce by s and c the sine and cosine 
solution of J (as a formal difference equation) , satisfying the initial conditions 

c(z,0) s(z,0) \ _ fl 0^ 
We observe that, we have that 

, . f c{z,n) s{z,n) \ A J_ ( ^ - ^ij) 1 

^•"•^^ 10(^^-1)^(^,^-1) a{n-l)s{z,n-l)^ 11 a(j) V «0')' « 

We note that 

, . det(z - J[o,„_i]) det(z - J[i,„_i]) 

(6.10) c(z, n) = „_i , s(^, ") = „n-i • 

11^=0 «U) 11^=0 au) 

For Im(z) > 0, we denote by u+(z, n) the solution of 

(6.11) Hu+^zu+, u+ef{Z+), u+(z,-l) = l. 
We then have that ii+(z,0) = — a(0)m+(z, J). We obtain that 

N 

(6.12) w+(z, iV) = (-1)^ Yi a{n)m+{z, /")). 

Tl = 

Hence, we obtain that 

Lemma 6.2. Assume (j6.3p . then for lm{z) > 

(6.13) I(^,{iV,}) = - hm -llog|7.+ (z,7V,)| 

1 



= - lim — logJ|ii+(z,iV,)|2 + |«+(z,iV, + l)|2 
= log(A^^) - / log\m+{z,{a,b))\dn{a,b). 



J 



Proof. Define L^{z, {Nj}) as the limit 

L+iz,{N,}) - - lim — J2 log|a(n)m+(z, J("))|, 
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which exists by continuity of J i— > log |to+(z, J)|, ()6.3|) . and (16. 7p . By (|6.12|) . 

\m+{z,J)\ < 

^\u+iz,N,)\^ + \u+{z,N, + 1)|2 = \u+iz,N,)\,Jl + HN, + l)mf ^+^^(z)|. 

we see that L+{z, {Nj}) is equal to the all the quantities on the right hand side of 
(|6.13p . In order to see the remaining inequality, observe that 

u+{z,N,) ^ TT J_fz~Kn) l\fm+iz,J) 

aiN, - l)u+(z, N, - 1)) a{n) { a{nf oj { 1 

which implies the claim by the Ruelle-Osceledec theorem. □ 

This shows (ii) of Theorem l6.ll Or next goal is to relate the Lyapunov exponent 
to the asymptotics of the cosine solution. Introduce for Im(z) > 

(6.14) mAr(z) = {Sn, {J[o,n] ~ 2;)"Mjv)- 
We have that 

Lemma 6.3. For Im(z) > 0, 

(6.15) ciz,N+l)_ 1 



c{z,N) a{N)mN{zy 
and \mN{z)\ < jK^- 
Proof. Introduce v for 0<n<iV+lby 

{0 n < 

{Sn,{Jlo,N]- 0<n<N 
^ n = N+l. 

One then checks that (J — z)v{n) — for 1 < n < and v{~l) — 0. Thus 
c{z,n) — C ■ v{n) for < 7i < iV + 1 for a fixed constant C. Now a computation 
shows that 

c{z,N + l) _v{N + 1) _ 1 
c{z,N) ~ v{N) ~ a{N)mN{zy 
which shows the claim. □ 

As in the proof of Lemma 16.21 one can now show that 



]- log |c(z, iV,) I = Ihn -1 log ^\c{z,N,)\^\+c{z,N^)\^ 



lim 

which implies for Im(z) > by the Ruelle-Osceledec theorem 

(6.16) lim ^\og\u+{z,N,)\ ^ - lim ^\og\c{z, N,)\, 

since the cosine solution c can never decay, since J is self-adjoint. In particular, 
the limit on the right hand side of (|6.16p exists for every z with Im(2:) > 0. 

Lemma 6.4. We have that 

1 f 1 

(6.17) -log|c(z,n)| = / log|z-t|dzy„- - Vlog|a(j)|. 

j=o 

Proof. This is a consequence of (16.101) and (|6.ip . □ 
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Lemma 6.5. Assume (|6.3|) . then 

(6.18) ly — lim vj^. 

j— »oo 

exists, and for Ini(z) > 

(6.19) lim -L log |c(z, N,)\ = log(A-i) + / log \t - z\dy. 
Furthermore, (j6.19p even holds for almost every z G R. 

Proof. (|6.18p follows from (|6.17p and the fact that the family of functions t 
log \t — z\ for Im(z) > separates points on the real axis. For the last statement, 
observe that (|6.17[) remains valid for z G K, and then use Theorem A. 7. in [13]. □ 

This shows (i) of Theorem 16. II Next, we observe that 
Lemma 6.6. For every E e R, we have that 

(6.20) L{E, {Nj}) < log(^-i) + j log \t - E\dv 

Proof. First observe that L{z, {-/Vj}) is a submean function of z, and z ^ log(A^^) + 
/log \ t ~ z\dv is subharmonic. This implies the claim by Theorem 1.1. in 0. □ 

We now come to 

Proof of Theorem \6.1\ (Hi). This is a consequence of the last lemma, and the fact 
that 



n 



1 



z — b{n) 



by (EH). 



□ 
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